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ABSTRACT 


A computer  data  reduction  system  has  been  reprogrammed  to 
sample  individual  scan  recordings  of  the  line- spread  function 
of  laser  light  propagating  through  the  turbulent  marine  boud- 
ary  layer. 

The  distribution  of  spatial  frequencies  of  the  Modulation 
Transfer  Functions  from  random  single  scan  samples  from 
experimental  data  sets  has  been  investigated.  A preliminary 
analysis  indicates  that  the  zero  spatial  frequency  component 
of  the  MTF  is  log-normal  distributed. 
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I.  INTRODUCTION 

Atmospheric  turbulence  adversely  affects  the  propagation 
of  a laser  beam  or  the  performance  of  an  imaging  system.  The 
Modulation  Transfer  Function  (MTF)  is  a measure  of  the  degra- 
dation of  resolution  and  thus  the  relative  ability  of  the 
atmosphere  to  transmit  spatial  frequency  information.  Knowl- 
edge of  the  MTF  of  the  atmosphere  is  necessary  for  prediction 
of  the  beam  cross-sectional  profile  on  target  and  spot  wander 
for  a high  energy  laser  and  its  beam  projection  system. 

Members  of  the  faculty  of  the  Naval  Postgraduate  School 
(NPS) , Department  of  Physics  and  Chemistry,  have  investigated 
the  effects  of  atmospheric  turbulence  on  the  propagation  of 

t 

optical  radiation  in  the  marine  boundary  layer. 

The  work  started  in  1973  and  has  continued  through  1978. 
The  research  has  been  concentrated  on  developing  a model  for 
the  MTF  depending  on  wavelength,  range  and  a properly  path- 
position  weighted  value  of  the  turbulence  structure  constant 
for  optical  index  Cn.  The  work  for  the  period  1 July  1974  - 
30  September  1977  is  described  in  Ref.  1.  Field  experiments 
were  performed  by  transmitting  laser  light  of  different  wave- 
lengths along  various  seapaths.  The  received  intensity  was 
detected  by  a slit- scanning  telescope  system  and  recorded  on 
tape.  The  Cn  value  was  determined  by  computer  calculations. 

The  objective  of  this  report  is  to  describe  the  develop- 
ment of  a computer  program  which  samples  individual  scan 
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recordings  of  the  line- spread  function  and  determines  a pre- 
liminary statistical  distribution  of  the  spatial  frequency 
content  of  the  Modulation  Transfer  Function. 


II.  ANALYSIS  OF  DATA 


A.  FUNCTIONAL  DESCRIPTION 

1.  Data 

Experimental  data  are  the  intensities  of  the  line- 
spread  functions  as  recorded  on  a 4-channel  magnetic  tape. 

One  track  is  used  for  verbal  identification  of  date,  time 
and  pertinent  experimental  information.  Two  tracks  are  used 
for  recording  scans  and  the  last  track  is  used  for  recording 
the  trigger  signal. 

Verbal  information  has  been  cross-checked  against 
experiment  log  book  and  MTF  data  output . Data  have  been 
recorded  for  one  and  two  wavelengths  simultaneously.  Only 
data  for  .6328  microns  have  been  analyzed  in  this  report. 

The  line- spread  function  was  sampled  by  the  scanning  detector 
at  a rate  of  50Hz  and  recorded  at  a tape  speed  of  37.5  inches 
per  second. 

2 . Data  Reduction 

The  data  reduction  is  done  on  a Nicolet  Instrument 
Company,  NIC-80  data  processor.  This  is  a data  digitizer 
followed  by  a 12K  20-bit  digital  computer  and  a Tektronix  CRT 
display  unit.  The  NIC-80  computer  is  controlled  by  a Hewlett- 
Packard  HP- 9825  desk  top  type  programmable  calculator. 

The  NIC-80  computer  reads  the  line-spread  function 
for  one  scan,  performs  a Fourier  Transform  (FT)  analysis  of 
the  signal  and  transfers  the  FT  spatial  frequency  components 
0 - 19  to  the  HP-9825.  The  MTF  of  the  atmosphere  is  calculated 
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by  dividing  the  FT  of  the  line- spread  function  by  the  FT  of 
the  optics.  The  FT  of  the  optics  is  stored  by  the  NIC-80  and 
transfered  to  HP-9825.  The  experimental  analog  signal  is 
digitized  at  1024  points.  The  DC  component  is  removed  by 
manual  setting  of  signal  bias  based  upon  observation  of  scan 
on  CRT  display.  Any  remaining  DC  component  is  compensated  by 
subtracting  the  average  intensity  value  of  32  points  at  each 
end  of  the  digitized  signal. 

When  not  performing  FT  and  transfer  operations,  the 
NIC-80  remains  in  signal  display  mode  with  successive  scans 
being  displayed  on  CRT. 

The  computation  and  transfer  time  for  one  scan  has 
been  observed  to  be  5 seconds,  and  will  be  the  time  between 
sampling  of  scans.  During  data  reduction,  the  tape  recorder 
was  operated  at  3.75  IPS,  corresponding  to  5 scans  per  second. 
This  implies  that  about  every  25th  scan  of  the  recorded  data 
was  sampled  for  analysis.  The  number  of  scans  sampled  was 
limited  by  the  length  of  experimental  data  recorded,  normally 
in  the  order  of  one  minute,  or  3000  scans.  This  allowed  a 
sample  size  of  100  - 125  for  data  analysis.  An  attempt  to 
run  the  recorder  at  .375  IPS  speed  was  not  successful  due  to 
problems  with  the  trigger  signal. 

3 . Statistical  Analysis 

The  HP-9825  calculates  the  MTF  at  20  different  spatial 
frequencies  from  0 to  38000  per  radian  in  steps  of  2000  per 
radian.  Thes**.  spatial  frequency  points  are  labeled  0 to  19 
and  are  tested  for  normal  and  log-normal  distribution.  Goodness 
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of  fit  is  determined  by  the  Chi  Square  test.  A data  summary 
containing  mean,  standard  deviation  and  Chi  Square  values 
for  normal  and  log-normal  tests  is  outputted  on  the  printer. 
Mean  and  + 1 standard  deviation  for  frequency  points  0-19 
is  plotted,  with  the  zero  spatial  frequency  normalized  to  a 
value  of  unity  and  all  other  data  adjusted  correspondingly. 

A histogram  showing  relative  distribution  of  MTF  values  for 
one  specific  spatial  frequency  is  plotted  as  requested  in 
program  input. 

Examples  of  the  above  are  given  in  Figures  1,  2 and  3 

respectively.  The  relationship  of  the  MTF  distribution  to 

the  Cn  value  has  been  investigated.  The  Cn  values  ranged  in 

-8  - 7 

the  order  of  8x10  to  1x10  . Typical  line-spread  functions 
for  two  different  Cn  values  are  shown  in  Figure  4.  Note  that 
Cn  squared  values  are  given. 

The  value  of  Cn  used  has  been  taken  from  existing  MTF 
data  outputs.  If  no  output  was  available,  an  MTF  analysis 
was  performed.  A detailed  discussion  of  MTF  and  Cn  is  given 
in  Ref.  1.  The  scan  presentation  is  a function  of  time  and 
detector  position.  The  dwell  time  setting  was  obtained  by 
analyzing  calibration  data  recorded  on  tape.  Dwell  time  is 
the  sampling  time  per  data  point.  Data  have  been  obtained 
by  recording  the  detector  output  with  a diffraction  grating 
of  known  dimensions  placed  in  front  of  the  telescope.  This 
provides  a known  reference  of  angular  spread.  The  dwell  time 
is  adjusted  so  that  the  angular  spread  of  the  first  order 
diffraction  as  observed  on  the  display  is  equal  to  the  ref- 
erence value. 
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relative  distribution  of 


UNITS  - 1/3  STANDARD  DEV  I RT 1 ON 

Figure  3 


Figure  4.  Typical  Line-Spread  Functions 


B.  PROGRAM  THEORETICAL  CONSIDERATIONS 

The  data  were  tested  for  normal  and  In  normal  distribution. 
1.  Normal  Curve 

The  normal  or  Gaussian  curve  is  expressed  as 
Y - constant  exp  (-%(^~ )^) 

S 

where 

x is  the  value  of  MTF  for  normal  test  distribution 
or  value  of  In  MTF  for  log-normal  distribution. 

x is  the  unbiased  estimate  of  the  Mean  of  the  MTF 
values  or  Mean  of  In  MTF  values . 

s is  the  unbiased  estimate  of  the  standard  deviation 
of  the  distribution. 

Y is  the  probability  density  of  the  distribution. 

Constant  = i/s/2iT,  where  i is  the  width  of  a histogram 
bin  to  make  area  under  curve  equal  to  i. 

The  Mean  x is  the  arithmetic  Mean  given  by: 

N 

x = E x. /N 
i-1  1 

where  N is  the  number  of  observations  (samples) . 

The  standard  deviation  s is  the  positive  root  of  the 

Variance  s squared  given  by 

o N _ 2 

s Z - E (x. -x)  / (N-l) 
i-1  1 

For  computational  purposes  the  following  algebraically 

equivalent  formula  has  been  used: 

, N 9 N 9 

» ( E (x.)z-(  E x. ) Z/N) / (N-l) 
i-1  1 i-1  1 
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The  probability  density  curve  has  been  computed  and 
plotted  over  a range  of  Mean  + 3 standard  deviations. 

2 . Goodness  of  Fit 

Test  of  goodness  of  fit  to  tested  distribution  has 
been  done  by  histogram  and  Chi  Square  calculation. 

The  width  of  the  histogram  bins  is  a program  variable 
expressed  as  a ratio  of  the  standard  deviation.  1/3  has  been 
used  corresponding  to  a total  of  18  bins.  The  number  of 
samples  with  values  within  a bin  is  expressed  as  a fraction 
of  1,  so  that  the  area  of  the  histogram  is  equal  to  the  width 
of  the  bin  interval. 

The  Chi  Square  goodness  of  fit  value  has  been  calcu- 
lated from 


where 

N = number  of  observations 
n ■ number  of  histogram  bins 
F(n)  = number  of  observations  within  a bin 

The  term  Nf(n)  is  the  theoretical  number  of  observations 
which  should  fall  within  a bin.  The  probability  area  has  been 
calculated  as  the  sum  of  a rectangle  and  a triangle  approxi- 
mating the  area  under  the  normal  curve.  The  sampling 
distribution  of  the  Chi  Square  statistics  is  tabulated  as  a 
function  of  degrees  of  freedom  and  percentile.  The  number  of 
degrees  of  freedom  is  calculated  from 
df  - k-n-1 
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where 

k - number  of  bins  in  histogram 

n - number  of  parameters  estimated 
For  a 18  bin  histogram  with  mean  and  standard  deviation 
estimated,  the  degrees  of  freedom  is  15,  with  a 95  percentile 
value  tabulated  as  25.  A calculated  Chi  Square  value  less 
than  25  will,  at  the  57.  significance  level,  imply  that  the 
data  fits  the  given  distribution. 

C . COMPUTER  PROGRAMS 

1 . NIC-80  Data  Reduction  Program 

This  program  has  been  developed  utilizing  subroutines 
from  the  existing  MTF  analysis  program.  Pages  3 and  4 of  the 
MTF  program  have  not  been  changed.  Main  program  is  on  page  1 
and  a data  transfer  subroutine  has  been  inserted  on  page  2. 

Program  flowchart  is  given  in  Figure  5A.  Program  listing  is 
on  page  122. 

2 . Data  Analysis  Program 

Due  to  limited  HP-9825  memory  capacity,  this  program 
has  been  split  into  two  parts  stored  on  separate  program 
files.  During  operation,  programs  are  loaded  and  executed 
under  program  control.  Program  "A"  is  on  file  8 and  contains 
NIC- 80  control  and  all  analysis  functions  except  data  summary, 
which  is  on  program  "B"  located  on  file  9.  Programs  are 
documented  on  flowcharts  (Figures  6 and  7) . Program  listings 
are  on  pages  117-121.  Execution  of  analysis  program  requires 
HP  tape  inserted  with  track  1 available  for  data  storage. 
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Each  data  file  must  be  minimum  160  bytes  long  and  the  number 
of  files  must  be  equal  to  or  larger  than  the  number  of  scans 
sampled. 

Execution  time  for  a 100  scan  analysis  run  is  of  the 
order  of  90  minutes.  Most  of  this  time  is  required  for  tape 
movements  and  operations . 


1 

f 

Calculate 
Mean  & 

S tandard 
Deviation 

1 

7 

Print 
MTF  0 & 
MTF  1 


Date:  6 figure  if,  date,  month,  year 

Time:  6 figure  if 

Cn:  As  argument  and  power 

Line  18 


For  MTF  spatial  frequencies  0-19 
Line  19-21 


To  verify  absolute  and  relative 

values 

Line  22 


Enter 
Class 
Interval 


Histogram  class  interval  is  expressed 
as  a ratio  of  standard  deviation,  e.g.  3 
Line  23 


or  already  LN  data  on  file 


Figure  6A.  Data  Analysis  Program 
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Calculates  and  plots  normal  curve 
based  upon  sample  estimate  of  mean 
and  standard  deviation 
Line  43-46 


Figure  6B.  Data  Analysis  Program 
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Ill . RESULTS 


A total  of  24  data  analysis  ms  have  been  performed. 

Five  of  these  have  been  on  the  same  data  sets  in  order  to 
investigate  consistency  of  distribution  within  the  same  data 
set. 

An  arbitrary  grading  system  for  goodness  of  fit  over  a 
range  of  spatial  frequencies  has  been  established  as  follows: 
All  components  in  range  within  95  percentile  ■ 1. 

One  component  in  range  outside  95  percentile  - 2. 

Two  or  more  components  in  range  outside  95  percentile 

» 3. 

Data  summary  and  MTF  mean  and  standard  deviation  plots 
are  given  on  pages  39-86,  arranged  in  order  of  decreasing 
order  of  Cn.  A distribution  analysis  summary  is  tabulated 
on  pages  31-33.  A set  of  histograms  for  all  spatial  frequen- 
cies for  one  data  set  is  given  on  pages  87-112- 

A.  CONSISTENCY  OF  DISTRIBUTION 

Data  set  number  12  has  been  run  5 times  to  verify  consis- 
tency of  distribution  as  only  every  25th  scan  was  sampled. 

Utilizing  the  grading  system  for  spatial  frequency  ranges 
0,  1-6,  7-12  and  13-19,  the  following  results  as  given  in 
Table  II  and  III,  were  found: 
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TABLE  I 


DISTIRUBTION  ANALYSIS  SUMMARY 


Set 

No. 

Cn 

Date 

Time 

ESHuol 

Norm 

Log 

Norm 

Comments 

1.04 

031176 

2320 

o 

1 

1 

0 freq.  log 

1 

1-6 

1 

1 

data  best  fit 

L 

h-» 

o 

1 

7-12 

1 

1 

13-19 

2 

3 

1.06 

021877 

1225 

0 

1 

1 

0 freq.  normal 

0 

1-6 

1 

1 

data  best  fit 

L 

10'7 

7-12 

2 

3 

13-19 

3 

3 

1.06 

040276 

1420 

0 

2 

2 

Large  standard 

3 

1-6 

3 

3 

deviation 

O 

1 

7-12 

3 

3 

13-19 

3 

3 

1.07 

031176 

2247 

0 

1 

1 

Good  fit  at 

4 

1-6 

1 

1 

low  frequen- 

t-* 

o 

1 

7-12 

1 

1 

cies 

13-19 

3 

3 

1.67 

031176 

2205 

0 

1 

1 

0 freq . log 

C 

1-6 

1 

2 

data  best  fit 

D 

o 

1 

'vj 

7-12 

1 

2 

13-19 

2 

3 

1.75 

040276 

1203 

0 

1 

1 

0 freq.  log 

6 

1-6 

3 

3 

data  best  fit 

10"7 

7-12 

3 

3 

13-19 

3 

3 

2.79 

092177 

1555 

0 

2 

H 

7 

1-6 

3 

S3 

oo 

• 

o 

t-» 

7-12 

3 

11 

13-19 

3 

■9 

4.25 

121776 

1140 

0 

2 

1 

log  data  best 

O 

1-6 

3 

1 

fit 

O 

i-1 

o 

• 

oo 

7-12 

3 

3 

13-19 

3 

3 

Legend:  I.  All  values  within  Chi  Square  95  percentile. 

2.  One  value  above  Chi  Square  95  percentile. 

3.  Two  or  more  values  above  Chi  Square 
95  percentile. 
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DISTRIBUTION  ANALYSIS  SUMMARY 


Set 

No. 

Cn 
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Log 

Norm 

Comments 
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1452 

IS 

1 

2 

0 freq.  log 

9 

00 
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O 
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2 

1 
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3 

3 
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3 

2 

5.62 

021877 
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0 

1 

1 

log  data  low 

10 

IQ'8 
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1 

1 

chi  sq.  values 

7-12 

2 

1 

13-19 

3 

2 

5.86 

121776 

1244 

0 

2 

2 

Absolute  values 

11 

00 

1 

o 

H 

1-6 

2 

1 

low 

7-12 

2 

3 

13-19 

3 

3 

5.89 

092377 

1522 

0 

2 

n 

1 

This  data  set 

12 

oo 

o 

i— i 

1-6 

3 

K3 

run  5 times  to 

7-12 

3 

check  consis- 

13-19 

3 

H 

tency  within 
data  set 

5.94 

092277 

1616 

0 

1 

2 

13 

1-6 

3 

3 

o 

■ 

oo 

7-12 

3 

2 

13-19 

3 

3 

092277 

1705 

0 

1 

1 

log  data  range 
0-8  low  chi  sq. 

14 

1-6 

1 

1 

7-12 

1 

2 

values 

13-19 

2 

2 

121776 

1012 

0 

1 

n 

log  data  good 

15 

1-6 

1 

is 

fit 

7-12 

3 

H 

13-19 

3 

a 

021577 

1700 

m 

2 

2 

16 

PI 

K59 

2 

3 

■29 

3 

3 

i VBtm 

11 

3 

Legend:  1.  All  values  within  Chi  Square  95  percentile. 

2.  One  value  above  Chi  Square  95  percentile. 

3.  Two  or  more  values  above  Chi  Square 
95  percentile. 
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DISTRIBUTION  ANALYSIS  SUMMARY 


Set 

No. 


17 


Cn 


Date 


Time 


Spat . 
Freq. 


Norm 


Log 

Norm 


Comments 


7.28 

io-8 


021777 


1900 


0 

1-6 

7-12 

13-19 


1 

1 

3 

3 


1 

1 

3 

2 


18 


19 


20 


7.35 

10"8 


021877 


1035 


0 

1-6 

7-12 

13-19 


1 

3 

3 

3 


2 

1 

3 

1 


Large  standard] 
deviations . 


8.23 

10'8 


021577 


1635 


0 

1-6 

7-12 

13-19 


2 

2 

3 

3 


2 

2 

3 

3 


Large  standard] 
deviations. 


8.56 

IO'8 


021577 


1600 


0 

1-6 

7-12 

13-19 


2 

3 

3 

3 


1 

1 

3 

3 


Legend:  1.  All  values  within  Chi  Square  95  percentile. 

2.  One  value  above  Chi  Square  95  percentile. 

3.  Two  or  more  values  above  Chi  Square 
95  percentile. 
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TABLE  II 


Grading  Consistency  of  Data  from  Same  Data  Set 


Frequency  Range 


Number  of  Runs  with  Grading 
Normal Log  Normal 


Grade 1 

0 1 

1-6  1 

7-12  0 

13-19  0 


The  number  of  samples  which  met  the 
for  each  spatial  frequency  in  the  0 
listed  in  Table  III. 


2 3 12  3 

4 0 4 1 0 

0 4 5 0 0 

0 5 3 1 1 

0 5 3 2 0 

95  percentile  criterion 
- 6 frequency  range  is 


TABLE  III 

Consistency  of  Individual  Frequency  Components 
from  Samples  of  Same  Data  Set 


Frequency 

0 

1 

2 

3 

4 

5 

6 


Number  of  Samples  within  957, 
Normal Log  Normal 


1 

1 

2 

3 

3 

2 

3 


4 

5 
5 
5 
5 
5 
5 
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I I 

The  results  indicate  that  the  distribution  can  be  expected  to 
give  same  Chi  Square  values  for  different  sample  series  from 
same  data  set. 

B.  FACTORS  AFFECTING  DISTRIBUTIONS 

1.  Scintillation 

If  the  intensity  of  a laser  beam  which  has  propagated 
through  a turbulent  medium  is  observed  at  a point,  the  inten- 
sity appears  to  fluctuate  or  scintillate  in  a random  way. 

It  has  been  shown  that  the  values  of  the  logarithm 
of  the  amplitude  follow  a Gaussian  distribution  (Ref.  1, 
page  15) . Spatial  frequency  component  zero  represents  the 
area  under  the  curve  (signal) . This  area  is  proportional  to 
the  received  intensity.  It  is  thus  expected  that  the  distri- 
' bution  of  the  zero  spatial  frequency  component  should  fit 
the  log  normal  distribution. 

2 . Signal  to  Noise  Ratio 

The  amplitude  of  the  spatial  frequencies  decreases 
with  increasing  frequency  number.  For  spatial  frequency  15 
and  above,  the  amplitude  is  of  the  order  of  107»  of  the  zero 
spatial  frequency  component.  At  higher  frequencies  noise 
will  contribute  more  relative  to  the  amplitude.  Noise  is 
considered  normal  distributed,  and  it  is  expected  that  the 
distribution  at  higher  frequencies  is  normal. 

C.  RESULTS  AND  DISCUSSION 

1 . Log-Normal  Distribution  of  0 Spatial  Frequency 

In  13  of  20  runs  the  spatial  frequency  zero  fits  the 


Ln  normal  Chi  Square  95  percentile  criterion.  For  runs  which 
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do  not  meet  the  distribution  criterion,  the  ratio  of  mean  to 
standard  deviation  is  in  the  order  of  2.  This  may  indicate 
saturation,  and  thus  not  follow  log  normal  distribution  (Ref. 
1,  page  16) . 

The  results  do  indicate  that  the  zero  spatial  freq- 
uency is  log  normal  distributed. 

2.  Normal  Distribution  of  Higher  Frequencies 

For  data  analysis  the  modulus,  e.g.,  absolute  value 
of  the  data  has  been  used,  and  sign  of  the  data  has  not  been 
included.  This  will  introduce  a skewness  of  distribution. 

It  is  believed  that  this  is  the  reason  for  lack  of  normal 
distribution  for  data  at  higher  frequencies. 

3 . Distribution  versus  Cn  Value 

Data  have  been  arbitrarily  grouped  into  ranges  with 
respect  to  Cn  values  as  follows: 

Cn  High:  Cn  > 6.10"8 

Cn  Medium:  Cn  between  2.10-8  and  6.10-8 

Cn  Low:  Cn  < 2.10-8 

The  ratio  of  distributions  graded  1 for  the  various 
Cn  ranges  is  listed  in  Table  IV. 

The  results  do  not  indicate  that  the  distribution 
does  change  with  Cn  values.  For  all  analysis  results,  it 
should  be  noted  that  the  sample  size  is  small,  and  the  results 
should  be  interpreted  as  indications  only. 
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TABLE  IV 


Ratio  of  Distribution  Graded 
I for  Different  Cn  Ranges 


Spatial 


Cn 

Frequency 

Normal 

Log  Normal 

H 

0 

5/6 

5/6 

I 

1 - 

6 

2/6 

3/6 

G 

7 - 

12 

3/6 

2/6 

H 

13  - 

19 

0/6 

0/6 

M 

0 

3/7 

4/7 

E 

D 

1 - 

6 

1/7 

5/7 

I 

U 

7 - 

12 

0/7 

3/7 

M 

13  - 

19 

0/7 

2/7 

0 

4/7 

4/7 

L 

1 - 

6 

3/7 

5/7 

0 

7 - 

12 

0/7 

1/7 

W 

13  - 

19 

1/7 

2/7 

D.  RECOMMENDATIONS  FOR  FURTHER  WORK 

To  reduce  analysis  computing  time,  the  HP-9825  should  be 
equipped  with  additional  memory  capacity.  To  eliminate  use 
of  tape  for  data  storage,  memory  should  be  increased  to  25K. 

The  data  analysis  for  each  run  should  be  stored  and 
analyzed  for  mean  and  standard  deviation.  Several  runs  on 
the  same  data  set  should  be  performed  to  verify  distribution. 
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Number  of  scans  sampled  in  each  run  should  be  increased 
by  either  solving  signal  trigger  problem  at  low  tape  speed 
or  by  reprogramming  NIC-80  computer  to  perform  and  store 
partial  analysis  results. 
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APPENDIX  A 

Operating  Instructions  for  Data  Analysis 
A.  Initialization 

1.  NIC-80,  Analysis  program  not  loaded 

(a)  Turn  on  computer  with  key,  depress  "STOP" 

(b)  Turn  on  Data  Bus  unit 

(c)  Set  toggle  switches  0 - 11  to  up  position 

(d)  Depress  "LOAD  PC" 

(e)  Depress  "CONTINUE" 

(f)  Press  "EXECUTE" 

2.  HP-9825 

(a)  Turn  on  HP-9825 

(b)  Insert  Data  Analysis  cassette  and  load  track  0, 
file  1 for  NIC-80  loading  program 

(c)  Push  "RUN";  c,  d,  i,  r,  n,  o,  1,  s,  t is  displayed. 

(d)  Push  "n"  and  "CONTINUE";  load  file  # is  displayed. 
The  NIC-80  data  analysis  program  is  located  on  files  2,  3,4 
and  5 . Enter  one  of  these  numbers  and  push  "CONTINUE" . 

(e)  Display  - c,  d,  i,  r,  n,  o,  1,  s,  t.  Push  "1" 
and  continue.  The  section  of  the  program  will  be  transferred 
to  NIC- 80.  At  the  end  of  each  transfer,  check  sums  are 
printed  and  should  be  identical. 

(f)  Repeat  steps  (d)  and  (e)  until  entire  program  is 
transferred.  If  NIC-80  contained  MTF  analysis  program  pre- 
viously, only  files  2 and  3 are  necessary  to  transfer. 
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(g)  On  NIC-80,  push  "STOP",  depress  "START"  and  push 
"EXECUTE"  on  completion  of  program  loading. 

3 . Tape  Recorder 

(a)  Turn  switch  to  3.75  IPS 

(b)  Switch  on  DC  bias  boxes . These  should  be  turned 
off  during  non-operating  periods  to  avoid  draining  batteries . 

(c)  Set  frequency  control  knob  on  appropriate  channels 
to  FM  and  lkc  bandpass  for  signal  and  trigger.  Voice  channel 
requires  37.5  IPS  and  100  Kc  for  operation. 

(d)  Turn  on  oscilloscope. 

4.  Pulse  Generator 

(a)  Turn  on  power 

(b)  Set  to  external  trigger  pulse. 

(c)  Set  width  to  10  microseconds.  On  tape  recorder, 
depress  "DRIVE"  and  adjust  oscilloscope  for  optimum  display. 

B.  EXECUTION 
1.  HP-9825 

(a)  Load  track  0,  file  8 and  push  "RUN" 

(b)  Data  entry; 

Wavelength,  e.g.,  .6328 
Optics,  e.g. , .4572 
Dwelltime,  e.g.,  95 
Number  of  scans,  e.g.,  125 

Push  "CONTINUE"  after  each  entry.  On  last  "CONTINUE"  NIC-80 
will  enter  Display/ Standby  mode.  HP- 9825  will  display 
"CONTINUE  WHEN  READY". 
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2.  NIC-80 

(a)  Verify  green  control  light  is  on 

(b)  Adjust  CRT  for  optimum  display  of  scans 

(c)  Adjust  DC  level  by  knobs  on  bias  boxes.  Use  low 
resolution  on  vertical  display  scale,  and  adjust  signal  to 
center  around  horizontal  scale. 

(d)  Adjust  horizontal  position  of  signal  by  pulse 
generator  delay  knob  so  as  to  center  around  vertical  axis . 

3.  To  start  data  collection,  push  "CONTINUE"  on  HP-9825. 
When  all  scans  are  completed,  NIC-80  will  return  to  Display/ 
Standby  mode  and  may  be  taken  off-line. 

4.  Data  Analysis  (HP-9825  only) 

(a)  On  completion  of  data  collection,  "HISTOGRAM  AND 
CONTINUE"  is  displayed.  Check  that  plotter  and  printer  are 
ready  before  pushing  "CONTINUE". 

(b)  Data  entry; 

Date,  e.g.,  021577  (month,  day,  year) 

Time,  e.g.,  1423 
Cn  arg , e.g.,  8.07 
Power,  e.g.,  -8 

Push  "CONTINUE"  after  each  entry. 

(c)  Unnormalized  values  for  MTF  spatial  frequencies  0 
and  1 are  printed  for  verification  of  absolute  and  relative 
values . 

(d)  Data  entry;  "No.  of  class  intervals"  is  displayed. 
Enter  number  of  histogram  bins  per  unit  standard  deviation. 
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_ ■ — — — 1 

Normally  3.  Push  "CONTINUE”. 

(e)  Display  "FREQ  PLOT,  YES-1".  Push  "1"  and  "CON- 
TINUE" if  plot  of  MTF  mean  + 1 standard  deivation  for  spatial 
frequency  0 - 19  is  required. 

(f)  Display  "LN  DATA?  YES-1".  Enter  1 and  push 
"CONTINUE"  if  data  analysis  of  normal  distribution  is  completed. 

Display  "LN  ANALYSIS".  Push  "CONTINUE"  if  data 
is  to  be  converted  to  Ln  data. 

(g)  Display  "DATA  SUM?  YES-1".  Enter  1 and  push 
"CONTINUE"  if  data  summary  is  required. 

Display  "SUMMARY,  CONT.”.  Push  "CONTINUE"  to 
start  execution  of  data  summary  program. 

(h)  Display  "SPAT  FREQ".  Enter  spatial  frequency 
number  (0  - 19)  for  which  histogram  will  be  plotted. 

Display  "REPLACE  PAPER  AND  CONTINUE".  Push 
"CONTINUE"  when  plotter  is  ready. 
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IBIS  PASS  IS  BEST  QUALITY  PRACTICABLE 
rods  copy  rural  shed  to  ddc  


0:  if  X-l;gto  "START" 

i:  mdec;par  l;fmt  l,fzl0.5;enp  "wavel" ,W, "oDt" ,0 
2:  enp  " DWELLTIME" , r 0 , "#  OP  SCANS", N 

3:  20*P;dim  F [P] , A [P] ,M (PJ ,S [P] ,X [P] ; wr t 2 .1, W, 0, r 0 , N; par  0 
4:  for  J-l  to  20;rdb(2)*A;rdb(2)*B;512*B+A*F[J] 

5:  if  F[J]>524287;F[J]-108576*F[J] 

6:  next  J;trk  1 

7:  enp  "CONTINUE  WHEN  READY", rO 
8:  wtb  2,r0 

9:  for  1*1  to  N;for  J*1  to  20 ;rdb(2)*A;rdb(2)*B;512*B+A*A[J] 
10:  if  A[J]>524237 ;A[J] -1048576*A[ J] 

11;  rdb  ( 2)  *A; rdb ( 2 ) *B ; 512*B+A*rO 
12:  if  r0>524287;r0-1048576*r0 
13:  A[J]*A[J]+rO*rO*rO;/rO*A[J] 

14:  A [J] /F  [J]*A [J] ;A[J]*A[J]+S[J]-S[J] ;M [ J ] +A [ J ] *M [ J] 

15:  next  J;rcf  I , A [ * ] 

16:  next  I 

17:  dsp  "HISTOGRAM  ?,and  continue" ; stp  ;flt  5 

18:  enp  "DATE" ,r 99 , "TIME" , r98 , "CN  ARG" , r94 , "POWER" , U 

19:  "PLOT"  :f or  Q*1  to  20;trk  1 

20:  (S[Q]-M[Q)*M[Q]/N)/(N-1)*S[Q]  ;M[Q]  /N-M  [Q  ] ; /S  [Q]  *S  [ Q] 

21:  next  Q 

22:  fit  5 ; pr  t "MTF0" ,M [ 1 ] , "MTF1 " ,M [2 ] ; if  r97*l;gto  "DS" 

23:  enp  "NO  CLASS  INTERV" , T;beeo 

24:  "START" : enp  "FREQ  PLOT, YES-1" , r9 2 ; if  r92*l;gto  "FREPLO" 
25:  if  r97*l ;g to  "DS" 

26:  "L" : enp  "LN  DATA?, YES*1 " , r97 ; i f r97* 1 ; 1*Z ;gto  "LN" 

27:  "DS":0»X;enp  "DATA  SUM?, YES-1" ,X ; if  X*l;trk  0;ldf  9 

28:  "PLOTO" :eno  "SPAT  FREQ" , A;A+1*Q; trk  1 

29:  dsp  "REPLACE  PLOT  PAPER  AND  CONTINUE" ; been ; stp 

30  : MlQ]-3*StQ]*C;M[Q]+3*S[Q)  -D 

31:  S[Q]/T*R;C-2.2*R*r91;csiz  1.6, 1.5,1 

32:  scl  C-2*R,D,-.2, 1.7; pit  D,0,l;plt  C-R/10, 0, 2 ;f xd  1 

33:  pit  r9  l+R/2 , 0, 1;  pit  r91  + R/2  , . 7 ,2  ;for  1-0  to  ,7  by  .1 

34:  ltr  r91+R/l.  5, 1;  lbl  I;plt  r91  + R/2-R/10, 1,1 

35:  pit  r91+R/2+R/10, I,2;next  I 

36:  -3*T*r  88; fxd  0 

37:  for  I*C  to  D by  R;olt  I,.01,l;plt  I, -.01, 2 

38:  itr  I-R/3 , -.08 ; lbl  r88 

39:  r88+l-»r88;next  I;pen;csiz  1.6, 1.5,1 

40:  ltr  C+3*R, 1 .6 ; lbl  "RELATIVE  DISTRIBUTION  OF  MTF" 


Figure  8A.  Data  Analysis  Program  Listing 
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41:  pen;plt  M[Q],0;plt  M[Q],.7,2 

42:  ltt  C+R,-.15;lbl  "UNITS  * 1/",T,"  STANDARD  DEVIATIONS" 

43:  "CURVE" :for  I*C  to  D by  (D-C) /100 ; (I-M [Q] )/S [Q]*r 1 

44  : rl*rl/2-rl;exp(-rl)*R/(S[QW(2*TT)  )t2 

45:  pit  I,r2,-2;next  I ;C*r 2; r2-R*r3 

46:  C*r2;r2+R*r 3;f It  5;A+A 

47:  for  K»ll  to  11+T*6 ; 0*rK; next  K 

48:  "HIST": for  I=»l  to  N;  C*r2 ; r2  + R*r  3 

49:  ldf  I,  A [*] ; for  K-ll  to  11+T*6 

50:  fit  5;  if  A[Q]>r2;if  A[Q]<r3  ;rK+l-rK 

51:  r3*r2;r3+R*r3;next  K;next  I 

52:  C*r2;r2+R*r 3; for  K»ll  to  11+T*6  ;rK-*r  1 

53:  rl/N*r l;pen;plt  r2,0,~2;plt  r2,rl,2;plt  r3,rl,2 

54:  pit  r3,0,2 

55:  (r  2-M  [Q ] )/S {Q] *r 7; r7*r 7/2+r7 ; exp  (-r7)  *R/(S  [Q]*/(2*ir))*r4 

56  : ( r 3— M [ Q ] )/S  [Q)*r7;  r7*r7/2-r7  ;exp  (-r7  ) *R/(S  [Q]*/  (2*ir) ) -r5 

57  : R*r4+abs  (r  5-r  4 ) /2*R-*-r 8 ; r8/R*r  8 

58 : (N*r8-rl*N)“2/(N*r8)-r9;r9+rlO-rlO 
59:  r3*r2;r3+R*r3;next  K 
60:  "LABEL" :C+R/3*r  1;  f xd  0;ltr  rl,1.45 
61:  lbl  "SPATIAL  FREQUENCY" , Q-l 

62:  if  r97*l; ltr  rl,1.25;lbl  "TEST  DIST  LN  NORMAL" ?gto  +2 
63:  ltr  rl, 1.25; lbl  "TEST  0I3T  NORMAL" 

64:  ltr  rl, 1.35; lbl  "NO  OF  SCANS", N 

65:  f xd  4 ; ltr  rl,1.15;lbl  "WAVELENGTH " ,W;  ltr  rl,1.05 

66:  lbl  "OPTICS", 0 

67:  fxd  0 ; ltr  rl,.9;lbl  " DATE " , r99 , " TIME" , r98; fit  5 
68:  M [Q] +R/3+r 1 ; ltr  rl,1.35;lbl  "St.  dev",S[Q] 

69:  ltr  rl,1.45;lbl  "MEAN" ,M [ Q] 

70:  ltr  r 1, 1 . 23;  lbl  "CHI  SQ",rlO;ltr  rl,1.15;fxd  3 
71:  lbl  "CN" ,r 94 ; f xd  0 
72:  lbl  "E" , U 

73:  for  I»1  to  90;0*rl;next  I 
74:  if  r97»l;qto  "DS" 

75:  gto  "L" 

76:  "LN":dsp  "LN  ANALYSIS" ; stp  ;trk  1 
77:  for  J*1  to  20;  G*M[  J)  *S  [ J)  ;next  J 
78:  for  I»1  to  N 

79:  ldf  I ,A  [* ) ; for  J-l  to  20;if  A[J]-0;if  J»20;gto  +5 
80:  if  A[J]»0;next  J 


Figure  8B.  Data  Analysis  Program  Listing 
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81:  ln(A[J]  )*A[J] 

82:  A[J]*A(J]+S[J]*S[J]  ;M  [ J]+A  [ J] -M  [J] 

83:  next  J;rcf  I,A[*] 

84:  next  I;gto  "PLOT" 

85:  "FREPLO" : scl  -3 , 20 ,-. 25 , 1 . 5 ; f xd  0;csiz  1.6, 1.5,1 
86:  pit  20 , 0 , 1; pi t 0,0,2;plt  -1.85 ,0 ,1 ;plt  -1.85,1,2 
87:  for  1*1  to  20;plt  1-1 , . 01 , 1 ;plt  1-1,- .01, 2 ; next  I 
38:  for  1*1  to  20  by  5 ; I-l»r95 ; 1 tr  r 95- . 3 , - .0 5; lbl  r95 
89:  next  I;fxd  1 

90:  for  1*0  to  1 by  ,l;ltr  -1.85,I;lbl  I;plt  -1,9, 1,1 
91:  for  1*0  to  1 by  ,l;ltr  -1.85,I;lbl  I;plt  -1.9, 1,1 
92:  pit  -1 .8 , I, 2 ; next  I 

93:  ltr  2, -.19; lbl  "SPATIAL  FREQUENCY  IN  1/  MILL I RAD IANS" 
94:  ltr  2,1.45; lbl  "MTF  MEAN  AND  + - 1 STANDARD  DEVIATION" 
95:  fxd  0 ;1 tr  2,1.35 

96:  lbl  "DATE"  ,r  99  , " TIME",r98,"  NO  OF  SCANS", N 

97:  fxd  4; ltr  5, 1.25; lbl  "WAVELENGTH", W, ” OPTICS ",0 
93:  fit  5 ;1 tr  5,1.15;fxd  3;lbl  "CN",r94;fxd  0;lbl  "E",U 
99:  M[l]*r93 

100:  csiz  3 , 1 , 1, 90 ; ltr  -2.1,.4;lbl  "MTF" 

101:  for  Q*1  to  20;50+Q-r48;70+Q-r49;M[Q]/r93-r43 
102:  S [Q] /r  93*r  49 

103:  pit  Q-l, r 48-r49 , 1; pit  Q-l , r48+r 49 , 2 

104:  pit  Q-l .05, r 43 , 1 ; pit  Q- . 9 5 , r4 3 , 2 ; pen ; nex t Q;gto  "L" 


Figxire  8C. 
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THIS  PACE  IS  BEST  QUALITY  PRACTICABLE 
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0:  dsp  "SUMMARY ,CONT" ;beeo ;stp 
1:  "X" : for  Q*1  to  20 
2:  S[Q]/T*R 

3:  M[Q]-3*S[Q]-C;M[Q]+3*S[Q]*D 
4:  for  K= II  to  11+T* 6 ;0*r K ; next  K 
5:  "HIST":for  1*1  to  N;  Or2  ;r2+R-*>r3 
6:  trk  1 

7:  ldf  I,Al*];for  K*ll  to  11+T*6 
8:  fit  5; if  A[Q)>r2;if  A [Q] < r 3 ; rK+l-r K 
9:  r3*r2;r3+R*r3;next  K;next  I 


10 

C*r  2 

; r 2+R+r 

3;for  X* 1 1 to  ll+T*6;rK* 

rl 

11 

rl/N 

*r  1 

12 

( c 2- 

M[Q] )/S 

[Q]*r7;r7*r7/2*r7;exD(-r 

7) 

*R/(S 

roj  * 

1(2* 

ff)  ) *l'4 

13 

(r  3- 

•MtQ])/S 

[Q]*r7; r7*r7/2*r7 ;exp (-r7 ) 

*R/(S 

[Q)  * 

V (2* 

n) ) *r5 

14 

R*r4+abs(r5 

-r4)/2*R--r8;  r8/R-r8 

15 

( N*r 

8-r 1*N) 

* 2/(N*r8) *r9; r9+rl0*rl0 

16 

r 3-*r 

2;  r3+R* 

r3;next  K ; rlQ+X  [Q)  ;Q-*r  10 

;next  Q 

;gto 

"PP 

INT" 

17 

"PRINT" :fmt 

2, 25x , "DATE ",2x,f6  ,0,2x 

1* 

9 

TIME" 

, 2x, 

f4 .0 

18 

fmt 

1, 2 5x , " 

Cn" , 2x,f 5 .3 ,2x, "E",f3 .0, 

2 x 

, "OPTICS" 

9 2 x # 

f 6 .4 , 2x 

19 

fmt 

3, 35x , " 

DATA  SUMMARY" 

20 

if  Z 

*l;gto 

"LI" 

21 

wr  t 

6 .3,13 

22 

wtb 

6,10,10 

23 

wr  t 

6 . 2 , r 9 9 

, r 98 

24 

wtb 

6,10,10 

25 

wr  t 

6 . 1 , r 94 

,U,0 

26 

wtb 

6,10,10 

27 

fmt 

4 , 2 5x , " 

TEST  DIST. NORMAL", 26x, "TEST  DIST  LN 

NORMAL" 

28 

wr  t 

6 .4 

29 

wtb 

6,10,10 

30 

fmt 

5, 15x, " 

SPAT  FR" , 5x , "MEAN  ", llx , " 

ST 

.DEV" 

, lOx 

, " CHI  SQ" , lOx 

31 

wr  t 

6 .5 ; wtb 

6, 10,10 

32 

fmt 

6 , 15x , f 

2.0,6x,ell.5,5x,ell.5,5x 

,e 

11.5 

33 

for 

I®1  to 

20 ; wr  t 6 .6,1-1, M(U  ,S[I) 

,x 

[ I]  ; next 

I 

34 

for 

1*1  to 

20; wtb  6,27, 10; next  I 

35 

fmt 

7 , 7 2x,e 

11.5 

36 

if  Z 

*1; gto 

"LI" 

37 

trk 

0;laf  8 

38 

"LI" 

:fmt  7, 

72x ,ell . 5 

39 

for 

1*1  to 

20;wrt  6 ,7  ,X  [ I ] ;next  I 

40 

trk 

0 ; ldf  8 

Figure  9. 
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0 

110600 

MEMA(6  00  /READOUT  LIGHTS 

1 

4304 

LCWORD  /CONTROL  WORD 

2 

44034 

STU  /READ  STATUS 

3 

10377 

ANDA ( 3 77  /COMPARE  SETTINGS 

4 

462340 

A— MZ (34 0 /CHECK  HARD  SETTING 

5 

0 

JMP  0 /RETURN  CHECK  LOOP 

6 

2175677 

ZERMA  1677  /0  TO  ACC  & 1677 

7 

4304 

LCW ( 0 /LOAD  CONTROL  WORD 

10 

3001737 

JMSGINP  /INPUT 

11 

3001745 

JMS01745  /PUTAC 

12 

1616 

1616  /ADR  OF  WAVELENGTH 

13 

3001737 

JMS3INP  /INPUT 

14 

3001745 

JMS91745  /PUTAC 

15 

1614 

1614  /ADR  OF  OPTICS 

16 

3001737 

JMS@ INP  /INPUT 

17 

3001745 

JMS91745  /PUTAC 

20 

1600 

1600  /ADR  OF  DWELLTIME 

21 

3001737 

JMSgTNP  /INPUT 

22 

3001745 

JMS91745  /PUTAC 

23 

1620 

1620  /ADR  NO  OF  SCANS 

24 

110020 

MEMA ( 20  /COUNT  VALUE 

25 

2405677 

AC CM  1677  /COUNT  INDEX 

26 

2001500 

JMS  1500  /CAL  ADRESS 

27 

3111576 

MEMA@ 1576  /CAL  VALUE 

30 

10777 

ANDA(777  /9  BITS 

31 

6474 

HSP  /IS  HP  9825  READY 

32 

31 

JMP  31  /WAIT  LOOP 

33 

4473 

PHS  /WRITE  TO  HP  9825 

34 

3111576 

MEMA?  1576  /CAL  VAL  RECALLED 

35 

405031 

RIS  11  /RIGHT  SHIFT 

36 

64  74 

HSP  /IS  HP  9825  READY 

37 

36 

JMP  36  /WAIT  LOOP 

40 

4473 

PHS  /WRITE  REMAINING  DATA 

41 

2125576 

MPOM  1576  /INCR.  CAL  ADRESS 

42 

2703677 

MMOZ  1677  /DECR.&TEST 

43 

27 

JMP  27  /RETURN  LOOP 

44 

3001742 

JMS31742  /GETAC 

45 

1600 

1600  /ADR.  OF  DWELLTIME 

46 

3001715 

JMS? 1715  /FIX 

Figure  10A.  Nic  80  Program  Listing 
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47 

3111722 

MEMA@1722  / FACM 

50 

4311 

LDWELL  /LOAD  DWELL  TIME 

51 

110000 

MEMA ( 0 /0  TO  ACC 

52 

4312 

LDELAY  /LOAD  DELAY  TIME 

53 

2111711 

MEMA  1711  /STAND  BY  LOOP 

54 

4304 

LCWORD  /LOAD  CONTROL  WORD 

55 

4302 

SET  /SET  MEASURE  FF 

56 

5220 

STOP  /READ  1 SCAN 

57 

6464 

HSRF  /SKIP  ON  HP  REQUEST 

60 

53 

JMP  53  /RETURN  STAND  BY  LOOP 

61 

44463 

RHSR  /READ  HP 

62 

462314 

A— MZ (314  /L  FOR  FUTURE  BRANCH 

63 

64 

JMP  64  /FUTURE  BR. ADR. 

64 

3001723 

JMS@ 172  3 /SHORT  TERM  RESET 

6 5 

2111704 

MEMA  1704  /B4 (0 ) 

66 

2405677 

ACCM  1677  / 

67 

2111626 

MEMA  1626  /#  OF  SCANS 

70 

2405672 

ACCM  1672  / 

71 

2111712 

MEMA  1712  /CONTROL  WORD  #2 

72 

4304 

LCWORD  /LOAD  CONT.WORD 

73 

4302 

SETM  /SET  MEASURE  FF 

74 

5220 

STOP  /READ  1 SCAN 

75 

100 

JMP  100  / 

76 

0 

0 /NOT  USED 

77 

0 

0 /NOT  USED 

100 

110200 

MEMA (200  / 

101 

4304 

LCWORD  /DISPLAY  CONTROL 

102 

2111700 

MEMA  1700  /30 ( 0 ) .LEFT  END 

103 

2405677 

ACCM  1677  / 

104 

2711701 

MMOA  1701  /B0 (1777) .RIGHT  END 

105 

2405676 

ACCM  1676  / 

106 

110040 

MEMA (40  /32  POINTS 

107 

2405675 

ACCM  1675  / 

110 

2165674 

ZERM  1674  / 

111 

3111677 

MEMA01677  /INTEGRAL  BASE  LINE 

112 

2505674 

A+MM  1674  / 

113 

3111676 

MEMA81676  / 

114 

2505674 

A+MM  1674  / 

Figure  10B.  Nic  80  Program  Listing 
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115 

2125677 

MPOM  1677  / 

116 

2705676 

MMOM  1676  / 

117 

2707675 

MMOMZ  1675  / 

120 

111 

JMP  111  / 

121 

2111674 

MEMA  1674  / 

122 

5026 

RASH6  /DIVIDE  BY  6 

123 

2405674 

AC CM  1674  /AVERAGE  DC  COMPONENT 

124 

2111700 

MEMA  1700  / B0 ( 0 ) 

125 

2405677 

AC CM  1677  / 

126 

3111677 

MEMA 016 77  / 

127 

2471674 

A-MA  1674  /SUBTRACTING 

130 

3405677 

ACCM01677  / 

131 

2135677 

MPOM A 1677  / 

132 

2463701 

A-MZ  1701  / 

133 

126 

JMP  126  /RETURN  ADR. LOOP 

134 

110002 

MEMA (2 

135 

3001722 

JMS01722  /MOVE  B0*B2 

136 

2111703 

MEMA  1703  /B3  ZERO  LOOP 

137 

2405677 

AC CM  1677  / 

140 

2111700 

MEMA  1700  /B0 

141 

2405676 

ACCM  1676  / 

142 

3165676 

Z ERM@ 16  76  / 

143 

3165677 

ZERM@1677  / 

144 

2125677 

MPOM  1677  / 

145 

2135676 

MPCMA  1676  / 

146 

2463701 

A-MZ  1701  / 

147 

142 

JMP  142  /RETURN  ADR. LOOP 

150 

3001776 

JMS01776  / 

151 

3001622 

JMS91622  /DPT  DATA  TRANSFER 

152 

4032 

ASCNTR  /ADVANCE  SWEEP  COUNTER 

153 

2703626 

MMOZ  1626  /DECREASE  SCAN  COUNT 

154 

71 

JMP  71  /RETURN  ADR. NEXT  SCAN 

155 

53 

JMP  53  /RETURN  ADR. STAND  BY  L. 

Figure  IOC.  Nic  80  Program  Listing 
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3445 

152 

152  /RETURN  ADRESS 

3446 

110024 

MEMA( 24 

/#  OF  SPAT.FR. 

3447 

2405675 

ACCM  1675 

/CONTING  INDEX 

3450 

2111702 

MEMA  1702 

/B2, REALS 

3451 

2405677 

ACCM  1677 

/CURRENT  B2 

3452 

2111703 

MEMA  1703 

/B3 , IMAGINARIES 

3453 

2405676 

ACCM  1676 

/CURRENT  B3 

3454 

3111677 

MEMA01677 

/VALUE  B2 

3455 

10777 

AND A (777 

/9  BITS 

3456 

64  74 

HSP  /HP 

READY 

3457 

14  56 

JMP  1456 

/WAIT  LOOP 

3460 

4473 

PHS  /WRITE  TO  HP 

3461 

3111677 

MEMA01677 

/RECALL  B2 

3462 

405031 

RIS  11  /RIGHT  SHIFT 

3463 

6474 

HSP  / 

3464 

1463 

JMP  1463 

/ 

3465 

4473 

PHS  / 

3466 

3111676 

MEMA@1676 

/ 

3467 

10777 

ANDA (777 

/ 

3470 

6474 

HSP  / 

3471 

1470 

JMP  1470 

/ 

3472 

4473 

PHS  / 

3473 

3111676 

MEMA@1676 

/ 

3474 

405031 

RIS  11  / 

3475 

6474 

HSP  / 

3476 

1475 

JMP  1475 

/ 

3477 

4473 

PHS  / 

3500 

2125677 

MPOM  1677 

/INCREASE  B2 

3501 

2125676 

MPOM  1676 

/INCREASE  B 3 

3502 

2707675 

MMCMZ  1675  /DECREASE  & TEST 

3503 

14  54 

1454  /RETURN  ADR. 

3504 

1001445 

JMP@ 1445 

/RETURN  TO  MAIN  PR 

Figure  10D.  Nic  80  Program  Listing 
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